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An action principle for the generalized harmonic formulation of general relativity is presented. 
The action is a functional of the spacetime metric and the gauge source vector. An action principle 
for the Z4 formulation of general relativity has been proposed recently by Bona, Bona-Casas and 
Palenzuela (BBP). The relationship between the generalized harmonic action and the BBP action 
is discussed in detail. 
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I. INTRODUCTION 

Einstein's equations can be expressed as an initial 
value problem using the familiar 3+1 splitting [H, 0|- For 
numerical applications, one must supplement the 3+1 
equations with coordinate conditions. Typically the full 
set of partial differential equations (PDE's) obtained in 
this way is not well posed (see, for example, Refs. @,ll|). 
Equations that are not well posed can be used for formal 
analyses, but they cannot be used for numerical appli- 
cations. Generalized harmonic (GH) gravity is a refor- 
mulation of Einstein's theory as a set of PDE's that is 
well posed. The GH equations are currently in use by a 
number of numerical relativity groups (see, for example, 
Refs. (^3). 

Einstein completed development of his general th eory 
of relativity in a series of papers published in 1915 |8l-[l0|. 
In the same year, Hilbert derived the field equations for 
general relativity by postulating a simple action principle 
motivated by general covariance [llj . The Hilbert action 
provides an economical and efficient way to define the 
theory. Throughout history, physicists have used varia- 
tional principles as a way of organizing and simplifying 
their descriptions of dynamical systems. Most physicists 
view the action as fundamental, and the classical equa- 
tions of motion as derived quantities. The action is typ- 
ically the starting point for a quantum analysis. 

Generalized harmonic gravity is a generalization of 
general relativity in the harmonic gauge. The general- 
ization to (in principle) arbitrary gauge conditions was 
first pointed out by Friedrich [12j . and later by Garfin- 
kle [13]. To my knowledge, the action for GH gravity 
has not been previously discussed. An action for general 
relativity in harmonic gauge was written down by Stone 
and Kuchar [3]. Their action was not complete in the 
sense that the harmonic coordinate conditions were not 
included among the equations of motion. Other efforts to 
write well-posed formulations of Einstein's equations in 
terms of a variational principle can be found in Refs. [1 5l — 
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Although generalized harmonic gravity is not a new 
theory, merely a reformulation of general relativity, the 
action principle presented in this paper provides a new 
perspective on the generalized harmonic system. This 
new perspective can help us understand the connection 
between GH gravity and other formulations of the Ein- 
stein equations. The GH action can serve as the basis 



for practical numerical calculations using variational or 
symplectic integrators [l8l - [20| . 

It is worth noting that any system of equations can 
be derived from a variational principle: Simply multi- 
ply each equation by an undetermined multiplier, add 
them together, and integrate over spacetime (for PDE's) 
or time (for ordinary differential equations). Such an ac- 
tion principle does not add any insights, and probably 
has no practical benefit. What we want in an action 
principle is an encoding of the equations of motion with- 
out the addition of extra unphysical variables that do not 
appear in the original differential equations. Not all sys- 
tems of equations can be derived from such a variational 
principle. For example, it appears that the Baumgarte- 
Shapiro-Shibata-Nakamura (BSSN) formulation of Ein- 
stein's theory [2l|, [22j cannot be derived from an action 
principle using only the BSSN variables. 

The action for GH gravity is presented in Sec. HU One 
of the features that emerges from this analysis is the need 
to introduce a background connection. The GH equa- 
tions are not usually written in terms of a background 
connection; equivalently, the background connection is 
usually set to zero. In numerical relativity applications 
this can be justified by choosing the background con- 
nection to be flat and interpreting the coordinates as 
Cartesian. Note that Kreiss, Reula, Sarbach and Wini- 
cour introduce a background metric in their studies of 
constraint-preserving boundary conditions for the gen- 
eralized harmonic equations 0, [24| ■ 

The Z4 system is a reformulation of Einstein's equa- 
tions that, with suitable coordinate conditions, is well- 
posed [HlHHl- Bona, Bona-Casas and Palenzuela (BBP) 
have recently proposed an action principle for Z4 p7| . 
In Sec. IIIII I discuss the relationship between the equa- 
tions of motion obtained from the BBP action and the 
Z4 equations, and point out their differences. The differ- 
ences are sublte and interesting. The key difference stems 
from the fact that the BBP action, like the familiar Pala- 
tini variational principle |27| , treats the spacetime metric 
and the connection as independent variables. As a result, 
the Ricci tensor that appears in the equations of motion 
for the BBP action is constructed from the independent 
connection and not from the Christoffel symbols. It is 
not clear whether or not the equations of motion for the 
BBP action have the nice properties of the Z4 equations. 
This shortcoming of the BBP variational principle can 
be corrected if we make a suitable change of variables 
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and eliminate the connection as an independent variable. 
The result is the GH action. 

In the Appendix I discuss the inverse problem of the 
calculus of variations. This provides a complementary 
perspective to the conclusions reached in Sec. IIIII In 
particular I argue that the equations of motion that fol- 
low from the BBP functional are not equivalent to the 
Z4 equations. A brief summary is contained in Sec. IIVI 



II. ACTION FOR GH GRAVITY 

Let g^v denote the spacetime metric and T a ^ denote 
the metric-compatible connection (the Christoffcl sym- 
bols). Let r Q Mi , denote a background connection that is 
torsion-free and therefore symmetric in its lower indices. 
We will use the shorthand notation 



Ar°, 



ra 
i 



ra 
i 



(i) 



for the difference between these connections. The sym- 
bol denotes the covariant derivative built from r CT M „. 
Note that Ar Q M „ is a type Q) tensor. Throughout this 
paper indices are raised and lowered with the metric <? M „ . 



Thus, for example, Ar 
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The generalized harmonic constraints are defined by 



C a = H a + Ar M / 



(2) 



where is the gauge source vector. The action for 
generalized harmonic gravity is the following functional 
of and 



(3) 



Here, is the Ricci tensor built from T 01 ^. Also, units 
have been chosen so that 16irG — 1, where G is Newton's 
constant. 



Before continuing, let me comment on the presence of 
the background connection. Since the Lagrangian must 
be a scalar density, then must be a covector. If we 
omit r CT MI/ from the definition ([2]), then must trans- 
form in such a way that + g ^ v g° '^T 1, ' a p is a covec- 
tor. Recall that under a change of spacetime coordinates, 
the transformation rule for the Christoffel symbols a p 
includes an inhomogeneous term. This inhomogeneous 
term, which is multiplied by g^g , must be canceled 
by a corresponding term from H^. It follows that the 
transformation rule for must include an inhomoge- 
neous term that depends on the metric. It is not possible 
for the transformation of to depend on the metric 
unless itself depends on the metric. However, for the 
moment, we would like to treat the metric g^ and the 
gauge source as independent variables in the action 
principle. For this reason, the background connection 
is needed to compensate for the inhomogeneity in the 
transformation rule for T a y, v . 

With the background connection included in the def- 
inition of the constraints C M , the gauge source is a 
covector. Although it is not logically necessary for 
to transform as a covector, as long as we are willing to 
give it a suitable dependence on g llVl it is at least con- 
venient for to transform as a covector. For example, 
we might find that a certain source works well for nu- 
merical simulations of black holes with a code that uses 
a Cartesian coordinate grid. Perhaps we would like to 
reproduce these results with a code that uses a spherical 
coordinate grid. If is a covector, we can easily deter- 
mine the correct form for the gauge source in spherical 
coordinates. 

Also observe that for most practical numerical appli- 
cations, it would be natural to choose T a M „ to be the flat 
connection. In this case the background connection com- 
ponents V 7 ^ u would be zero in Cartesian coordinates, but 
nonzero in spherical coordinates. 

Now consider the variation of the action ((3]). The func- 
tional derivatives of S[g^, H^] are 



SS 
SH^ 

SS 
5g^u 



-gc>< , 



(4a) 
(4b) 



where G^ v = R^-Rg^/2 is the Einstein tensor. Paren- 
theses around indices denote symmetrization. Note that 
V M is the covariant derivative built from the Christof- 



The background connection Y a ^ v appears in the action as an 
external field and is not varied. 



fel symbols r a Mt/ . It is related to the background co- 
variant derivative by V M V^ = V M K — Ar'^Vv, which 
holds for any covector V^. The vacuum Einstein equa- 
tions are obtained by setting the functional derivatives 
(|4]) to zero. Equation (]4^) tells us that C p = 0; hence 
C are constraints for the generalized harmonic system. 
With C' 1 =0, Eq. (|3Jd) reduces to the vacuum Einstein 
equations G^ v = 0. Matter fields can be included in a 
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straightforward way. 

A convenient form of the equations of motion is ob- 
tained by choosing \J—gg^ v and — y/— gH^ as indepen- 
dent variables in the variational principle, rather than 



6V 



and H^. This leads to the vacuum equations 
6S 



= 



5S 



S{V=gg' a ') 



-w) - C »> 



(5a) 
(5b) 



Note that the generalized harmonic equation are usu- 



ally written in the form R 



0. Neither 



Eq. ijlj)) nor Eq. (JSJd) is identical to the usual equation. 
The differences are terms proportional to the constraints 
C M . These terms depend on the choice of independent 
variables and are not particularly important. As we will 
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see, the presence or absence of these terms does not af- 
fect the properties that makes the generalized harmonic 
equations useful. 



The equations of motion ((4]) are equivalent to Ein- 
stein's equations. Of course, this assumes that each equa- 
tion holds for all time. In particular, the constraints 
C M = must hold for all time. We would like to re- 
interpret these equations as an initial value problem. For 
this purpose we follow the analysis of Lindblom, Scheel, 
Kidder, Owen and Rinne and derive two key results 



from Eq. (0Jd). Let n M denote the unit normal to a fo- 
liation of spacetime by spacelike hypersurfaces, and let 
hfiu = g^ + n^n,^ denote the metric induced on these hy- 
persurfaces. The first result is obtained by contracting 
Eq. ((3Jd) with n v , which yields 



G» v n„ - Vv CT C^ = \(W°n p - h a p n»)V a C p 



C&AT^Pg - C a AY^ v - -CT + -cTCX" 
h 2 4 



(6) 



The second result is obtained by letting the covariant derivative V u act on Eq. ^p) and using the Ricci identity. This 
gives 



V CT V CT C' 



(7) 



where the term V^G* 4 " has been set to zero by the con- 
tracted Bianchi identity. 

The first term on the left-hand side of Eq. ^ is the 
Hamiltonian and momentum constraints, which we de- 
note M.^ = G^n u . The second term on the left-hand 
side is proportional to n a \7 a C^ = {d t C^ - (3 l diC^)/a + 
n°T M „C '. Each of the terms on the right-hand side 
of Eq. © is proportional to the constraints C p or their 
spatial derivatives. It follows that Eq. ([6]) has the form 

M" - = {terms ~ C, d.C) , (8) 

where diC denotes spatial derivatives of C p . 

Now consider the initial value problem. Equation ([8jl 
tells us that if C M and Af vanish initially, then 9(C M van- 
ishes initially. Then Eq. (JT)) implies that C 1 will remain 
zero throughout the evolution defined by Eq. ((Dd). In 
turn, Eq. © tells us that AP will remain zero through- 
out the evolution. The same conclusion can be reached 
by splitting the derivatives in Eq. (JT)) into space and time. 
Together with Eq. ([5]) one finds the results 

d t C = {terms - M, C, d t C) , (9a) 
9 t M" = {terms ~ Al, d,M, C, d l C, d l dfi} . (9b) 



These equations are consequences of Eq. (J4Jd) alone. 
Therefore, if the constraints C and AP vanish initially, 
then the evolution equation (|3Jd) will maintain the values 
C M = AP = throughout the evolution. 

Observe that Eqs. (J4]b) and (|5Jd) are not equivalent. If 
we take the trace reversed version of Eq. (JSJj) and raise its 
indices, the result differs from Eq. ^p) by terms that are 
linear and quadratic in the constraints C M . The difference 
does not depend on derivatives of the C's. As a result, 
the arguments that led to Eqs. ^ hold for the evolution 
equation ^jp) as well. In fact, we are free to drop any 
terms in Eqs. (UJa) or (|5Jd) that are linear or quadratic in 
the constraints. 

The discussion above shows that the relations ^ hold 
for any equation of the form 



R»» - V ( „C U) = {terms ~ C) . (10) 



The terms proportional to C M can include, for example, 
constraint damping terms. From the definition of the 
Ricci tensor we have 



I 

+ g a ' 3 [-Ar^Arv + 2Ar CT Q(p Ar 1/)(3a + atvaiVwj] , (n) 
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where R° is the Riemann tensor built from the background connection T a p U . Then the evolution equation (p~0|) 
becomes 

+2g al3 [-AT aaP AT% lu + 2AT° a{f ,AT„ )l3rT + AT\ a AT a ^] + {terms ~ C) . (12) 



This is a wave equation for each component of the space- 
time metric. The initial value problem for the GH system 
is described as follows: Specify initial data for g^ and Hp 
that satisfies = A'P = 0, then evolve the metric with 
the wave equation (fT2"[) . Observe that the gauge source 
vector Hp is freely specifiable, apart from the restriction 
C M = at the initial time. 



we first solve Eq. ([TBI :) for T a p U . This equation can be 
rearranged to give 

2g^Z a = Q," p „g^ - 2^ v \ + 5^^% . (17) 

By setting v — a we obtain 



(18) 



III. BBP ACTION 



The functional 

S\g ia ,,Z ll) T a ia ,] = Jd i xV^99 tl " [iV + 2V M Z„] (13) 

was proposed by Bona, Bona-Casas and Palenzuela 
(BBP) in Ref. [ljj as an action principle for the Z4 for- 
mulation of general relativity. This action is a functional 
of the spacetime metric g^, a covariant vector Zp, and 
a torsion-free connection T a p U . The covariant deriva- 
tive V M is built from this connection. Likewise the Ricci 
tensor that appears in the Lagrangian is defined by 

R^v ^cjT p V c^r p G + r^^^r p a r^oT V p . (14) 



(This definition differs slightly from that of Ref. [17| . As 
defined here, Rp V is not necessarily symmetric.) We will 
frequently use the abbreviation 



QCT _ pCT _ p(7 

^ L flu — x fiv 1 fiv : 

for the difference between the connection T a 
Christoffcl symbols T a ^ 



(15) 
and the 



Note that indices are raised 
and lowered with gp V and its inverse. Thus, for example, 
f = a f p 

Variation of the BBP action yields the vacuum equa- 
tions 

SS 



= 
= 
= 



s(V^gg^) 
1 ss 



^SZp 

1 ss 



-g 5T° 



= W pa g» v - 2VL { ^ a 
2Z a g» v 



per 1 



(16a) 
(16b) 

(16c) 



For convenience, we have chosen the independent vari- 
ables to be \J—gg^ v , Z fll and P 7 ^,,. It is straightforward 
to include matter fields. 

Equations ([T5)l . which follow from the BBP action, are 
equivalent to Einstein's general relativity. To show this, 



where the spacetime dimension is assumed to be 4. Now 
take the trace over the indices \x and v in Eq. (fTT| to 
yield 

«V - . (19) 

Putting the results (IT7l - [T9|) together gives 

£V<r + &vn<7 = g (Z a g^ + Z{n9v)<r) ■ (20) 

Now write down two more copies of this equation with 
index replacements fi v, v a, a fi in the first 
copy and /i — >• a, v — >• /i, a — >• v in the second. Add 
the second copy to Eq. (|20l) , then subtract the first copy. 
This yields 



for the solution of Eq. ([16b ). 

The vacuum equation of motion (|16b ) implies 



Wp = 



(21) 



(22) 



With the result ([21]) we see that Eqs. ([TBp) and ([16b), 
together, have the solution 



Zp = 

nv = 



(23a) 
(23b) 



The second of these equations tells us that the connection 
Y a p U is equal to the Christoffel symbols. The results ([2"3"]l 
show that the equation of motion (ITBk ) is equivalent to 
the vacuum Einstein equations, Rp V = 0. 

The Z4 equations are usually written as Rp V + 
2V (A1 Z„) = and Zp = 0. The equation Rp V +2V( M Z 1/) = 
has the same key properties as Eq. ^3Jd) or ([5b>) for the 
GH system. By the same analysis that led to Eqs. ([9]), 
one can show that the equation Rp V + 2V( AI Z„) = im- 
plies 

dtZ^- = {terms - M, Z, d t Z} , (24a) 
d t M" = {terms ~ M, diM, Z, diZ, d^Z} .(24b) 



5 



Thus, if Z^ = and M.^ = initially, then Z M and 
.M^ will remain zero throughout the evolution defined 
by + 2V (Al Z„) = 0. 

Unfortunately, the equation (ITBb ) that comes from 
the BBP action does not appear to have this property, 
for two reasons. First, the trace-reversed Ricci tensor 
G^v = R^v — g^ v R a j3g a ^ /2, built with the connection 
r%„, does not satisfy the contracted Bianchi identities. 
Second, the Hamiltonian and momentum constraints are 
not equivalent to the normal projections of G^ v . The 
argument showing that Z^ and will remain zero, as- 
suming they are zero initially, does not obviously hold 
for the equation -ff( p „) + 2V( AI Z„) = 0. 

Since Eqs. (fTBb ) and (|T6b ). together, imply Y a ^ = 
T a M „ and Z M = 0, we are free to replace the connection 
f with the Christoffel symbols T a ^ when solving the 
equations of motion. It follows that Eqs. ([To]) are equiv- 
alent to the system 



-2Z a g 



I'" 



(25a) 
(25b) 



obtained by setting f = T a ^ in Eqs. (fTHk ) and (fTBb ). 
These are the Z4 equations. However, these equations do 
not appear to coincide with the extrema of any action 
functional. In other words, there is no functional (to 
my knowledge) whose functional derivatives are linear 
combinations of i? M ^ + 2V( A ,Z ;y ) and —2Z a g^ v . This point 
is discussed more thoroughly in the Appendix. 

Note in particular that the functional obtained by set- 
ting f "pi, — r CT M k in the BBP action (fl~3"|) docs not yield 
Eqs. (|25|) for its extrema. This is an example of a gen- 
eral rule: One cannot always reduce an action principle 
by using results from the equations of motion. Consider 
an action S[u, v] that is a functional of two sets of vari- 
ables, u % and v a . If the equations of motion 8S/8u l = 
can be solved for the variables u 1 as functions of v a , then 
it is indeed permissible to use the solutions u l — u l (v) 
to eliminate u l from the action. On the other hand, one 
or more of the equations 8S/5u l = might yield, for ex- 
ample, v 1 as a function of the other v's and the it's. It 
is not permissible to use this result to eliminate v 1 from 
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the action. 

In light of these remarks, consider the BBP action (Ti~3"]) 
and the equations of motion ([16]) . As the result (f2"Tj) 
shows, the equation (|16b) has the solution 



rcr -pa 



flu 



4 



(26) 



Ofor 



In this case we have solved the equation 5S/ ST a i 
T a M „ and we are allowed to use this solution to simplify 
the action. The result is 



Rav 



-Z,,Z,. 



and the equations of motion become 







SS 



1 5S 8 ,. 
\J—g oZu 3 



(27) 

(28a) 
(28b) 



These equations are physically correct — they are equiva- 
lent to vacuum general relativity. They do not, however, 
have the form of the usual Z4 equations. 

Another option is to solve the equations of motion 
flUjb) and ([TBb). together, for f CT M „ and Z„. The so- 



lution is listed in Eqs. 



If we use these results to 



eliminate T a ^ and Z^ from the action we are left with 
the Hilbcrt action. The equations of motion are the vac- 
uum Einstein equations which are, of course, physically 
correct. However, they are not the usual Z4 equations. 

The equation of motion (|16b ). by itself, does not imply 
V" fj, u — r CT due to the presence of the fields Z^. We 
can try to eliminate Z^ from the functional derivative 
8S/8T a ^ by changing the independent variables in the 
action principle. Since a change of independent variables 
will merely mix the equations of motion, it will not be 
possible to eliminate Z^ from SS/8T a ^ unless Z^ ap- 
pears undifferentiated in one of the other equations of 
motion. With a simple modification of the action, the 
fields Z^ will appear in the functional derivatives SS/SZ^. 
Thus, let 



S[g^v, Zfj,, r fiv] — 

so that the equations of motion become 

5S 



= 
= 
= 



1 ss 



LV \Rp,v + 2V M Z y + XZ^Z^ 

R(tiv) + 2V ' ( M Z v) + xz^z u , 
-2^% + 2\Z P , 

: Q p pa g» v - 2fi^ CT + 5^n v)p p - 2Z a g^ v 

I 



(29) 

(30a) 
(30b) 
(30c) 



Here, A is a constant parameter. 



We can now look for a change of independent vari- 
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ables that will mix the equation of motion (|30b) with 
P0b) , and in the process eliminate from the functional 
derivatives SSjSY" p V . This_is accomplished by replacing 
with a combination of T a ^ and a new independent 
variable, a covariant vector that we call Hp. For example, 
we can replace Zp with the linear combination 



(31) 

in the action (|29| . The resulting equations of motion are 
SS 



= 
= 
= 



1 SS 



R { p v) + \ZpZ v + {terms ~ JF^a) 



^gSHp 
1 SS 



-g arv 



= --W y a + 2Z» , 
A 



-O p r/ Il/ 



(32b) 



(32c) 



with Z p given by Eq. (|3T|). _Eq. ([3"2b) has the de- 
sired property — its solution is T a p V = T a p V (assuming 
A 7^ 4/3). However, Eq. ([32k ) no longer includes the term 
proportional to V^iv) that characterizes the Z4 equa- 
tion (125k ') . This is because the change of variables (|3"Tj) 
contains derivatives of the metric through the Christoffcl 
symbols. 

We can eliminate the Christoffel symbols V p V from 
the change of variables f|3 1 1) by replacing them with a 
background connection T a t2l/ . Therefore, let 



(ffp+fVViyv) (33) 



in the action (|29[) . The equations of motion become 





= 
= 



SS 



-99 
SS 



,,,'1 



-5 ^ 



2Z^ 



+ {terms - (ft^ - XZP)} 



-9 ST- 



(34a) 
(34b) 
(34c) 



r 



where Zp is given by Eq. (|33|) . The solution of Eq. (|34b ) 
is f? ^ = for A ^ 4/3, and we are allowed to use 
Y" a V = ^ a fj,u in the action to eliminate r CT ^„. In the 



process, the definition (f3"3"f becomes Z M = C M /A, where 
Cp is the generalized harmonic constraint ([2]) . The action 
becomes 



XV ~99 



**>LLU 



(35) 



When A = — 2 these are the GH equations and Eq. (|35j) 
is the GH action. 

The preceding analysis shows that we are naturally 
led to the GH action when we attempt to reformulate 
the BBP action without the connection T" ^y. The GH 
action Q can be obtained directly from the BBP action 
(fT3"l) by the change of variables 



where the term proportional to ^ pZ v has been integrated 
to the boundary and discarded. The equations of motion 
are 



SS 2 



i 



-C c 



= 



+ {terms ~ C CT } , 



1 SS 



(36a) 
(36b) 



(37) 



With this definition, the BBP action becomes 



S[giMv,H^,,T p V ] 



d A x- 



-99" 



(38) 



The equation of motion SS/SV 7 ^ = has the solution Substituting this result into the action (|38p and discard- 



(39) 
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ing a boundary term yields the GH action (|3]). 



IV. SUMMARY 

The action for the generalized harmonic formulation 
of general relativity has the remarkably simple form dis- 
played in Eq. (|3|) . This action can be used as the starting 
point for further developments, such as the Hamiltonian 
formulation of GH gravity. We can also use the action to 
develop variational and symplectic integration schemes. 
The BBP action presented in Ref. [ijj is closely related 
to the GH action, but the equations of motion that follow 
from the BBP action are not obviously equivalent to the 
Z4 equations. After a change of variables, the indepen- 
dent connection T" ^ can be eliminated from the BBP 
action, reducing it to the GH action. 
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Appendix: The inverse problem of the calculus of 
variations 

The problem of finding an action for the GH (or Z4) 
equations is an example of the inverse problem of the cal- 
culus of variations. This subject has a long history [28[. 
In its most basic form, the inverse problem of the calcu- 
lus of variations can be stated as follows. Given a set of 
differential equations E A (<fi, d<j>, . . .) = for the variables 
4> A , does there exist a functional S[<f>] whose functional 
derivatives are E A (cf), d<j), . . .)? If so, is the functional 
unique? The index A runs from 1 to N and d(j) represents 
the partial derivatives of the dependent variables (f> A with 
respect to the independent variables. For ordinary differ- 
ential equations, there is only one independent variable; 
for partial differential equations, there are two or more 
independent variables. The dots in E A (cf), dcf>, . . .) repre- 
sent higher order derivatives of <f) A . 

An acceptable action functional for the GH or Z4 equa- 
tions does not need to reproduce the differential equa- 
tions identically. It is acceptable if the functional deriva- 
tives of the action are a linear combination of E A . This 
formulation of the inverse problem of the calculus of vari- 
ations is often referred to as the variational multiplier 
problem [H [Hj]. Thus, given a system E A ((/), d<j>,...) = 
0, we seek a functional S[4>] that satisfies 



where M AB is an invertible matrix that depends on <j> A 
and its derivatives. Equation (IA. 1[) says that the ex- 
pressions E A are linear combinations of the functional 
derivatives of S [</>]. 

The inverse problem of the calculus of variations as- 
sumes that the action is a functional only of those vari- 
ables 4> A that appear in the system of equations E A = 0. 
(It also assumes that the number of equations is equal to 
the number of variables.) As an alternative, consider the 
functional S[cj),A] = / A A E A (cj), <9</>, . . .) of <j) A and A A . 
The functional derivatives of S[4>, A] include E A . Equiv- 
alently, the conditions for the extremization of S[(f>, A] 
imply E A — 0. In spite of this fact, the functional S[<p, A] 
is not considered a valid action for the equations E A = 
because it depends on the extra unphysical variables Aa- 

In the variational multiplier problem (|A. 1|) . M AB can 
depend on the fields (f> A and their derivatives but it is 
not allowed to be a differential operator. This restriction 
on M AB is a natural one, since we want the functional 
derivatives of the action to yield the same system of dif- 
ferential equations as defined by E A = 0. A derivative 
operator in M AB can change the differential order of the 
functional derivatives so that the extremum of the action 
is no longer equivalent to the original differential system. 
Although this can happen when M AB contains differen- 
tial operators, it does not always happen. 

Let us consider the consequences of this restriction in 
the context of the BBP functional (fT3|). The functional 
derivatives of the BBP action are displayed in Eqs. (|16[) . 
A close examination of the analysis following these equa- 
tions shows that the functional derivatives (flBb ) and 
(|16b ) can be rearranged, by a linear transformation, to 
form the left-hand sides of Eqs. ([23| . In other words, 
there is a matrix M AB that mixes the functional deriva- 
tives of the BBP functional, leading to the result (using 
matrix notation in place of the indices A and B) 



Mi 



5S 



(A.2) 



We can use the definitions (114]) and (115)) to write this 
result in the form 



Mi 



5S 



r^ 7 - r^ 7 



(A.3) 



For simplicity, only two terms are displayed in the first 
row. 

Now we ask whether there exists a further mixing of 
the functional derivatives that will yield the Z4 equations 
R^u + 2V( Al Z !y ) = and Z^ = 0. The mixture must 
replace derivatives of the background connection T a 
with derivatives of the Christoffcl symbols r CT „„ in the 
first row of Eq. (|A.3|) . The matrix that does this is 



M AB {$, 



SSI 



E' 



(A.l) 



i o ^)0a-az*SA) + -- 

Mi = | 1 
1 



(A.4) 
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where each of the l's is an identity tensor. In this ex- 
ample both Mi and M 2 are invertible. But because M 2 
contains a derivative operator, the matrix MiM\ does 
not qualify as a valid variational multiplier for the in- 
verse problem of the calculus of variations. The conclu- 
sion is that the BBP functional (p~3|) docs not qualify as 
an action principle for the Z4 equations. 

In the present example the differential operator M^M\ 
is invertible, and it does not change the differential order 
of the functional derivatives of the BBP action. So per- 
haps the restriction that M AB should not contain any 
derivative operators is too severe. Perhaps the only re- 
striction on M AB should be invertibility. Note, however, 
that if we allow M AB to be a differential operator then 
there exist action functionals for the Z4 equations that 
are more simple than the BBP functional. For example, 
the action of Eq. (|27| has functional derivatives 



with the invertible matrix 



6S 



(A.5) 



as seen from Eqs. (|28|) . These can be rearranged to give 
'SS\ ( + 2V (M Z I/) 



M 



S({> 



(A.6) 



M 



1 -(2 5 ^^)+3 5(7(M V 1/) )/(4 A 
-3W(8V=5) 



(A.7) 



If we allow M AB to mix SS/5cj) A with derivatives of 
SS/Scj) A , then by this criterion the functional (pZ7|) would 
be a valid action principle for Z4. 

The view among researchers who study the inverse 
problem of the calculus of variations is that the varia- 
tional multiplier should be an invertible matrix that de- 
pends only on the variables and their derivatives [28| . 
According to this view, neither the BBP functional ([13")) 
nor the functional of Eq. (|27[) qualify as action principles 
for Z4. The GH functional ©, on the other hand, is a 
valid action principle for the GH formulation of general 
relativity. In particular, the equations of motion Q or 
(O follow directly from this action and have the desired 
properties discussed in Sec. [TTJ 
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